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ABSTRACT 

We examine the velocity width of cool X-ray emitting material using XMM-Newton Reflec- 
tion Grating Spectrometer (RGS) spectra of a sample of clusters and group of galaxies and 
elliptical galaxies. Improving on our previous analyses, we apply a spectral model which ac- 
counts for broadening due to the spatial extent of the source. With both conventional and 
Markov Chain Monte Carlo approaches we obtain limits, or in a few cases measurements, of 
the velocity broadening of the coolest X-ray material. In our sample, we include new obser- 
vations targeting objects with compact, bright, line-rich cores. One of these, MACS J2229.7- 
2755, gives a velocity limit of 280 km s -1 at the 90 per cent confidence level. Other systems 
with limits close to 300 kms" 1 include A 1835, NGC4261 and NGC4472. For more than 
a third of the targets we find limits better than 500 kms" 1 . HCG62, NGC 1399 and A3112 
show evidence for ~ 400kms~ 1 velocity broadening. For a smaller sample of objects, we 
use continuum-subtracted emission line surface brightness profiles to account for the spatial 
broadening. Although there are significant systematic errors associated with the technique 
(~ 15010ns- 1 ), we find broadening at the level of 280 to 500 kms" 1 in A3112, NGC 1399 
and NGC 4636. 
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1 INTRODUCTION 

Measuring the velocity distribution of the intracluster medium 
(ICM) in clusters and groups of galaxies and in elliptical galaxies 
is important for understanding how such structures form and the 
strength of active galactic nucleus (AGN) feedback. Simulations of 
AGN feedback i mply the creation of ICM velocity widths of several 



hundred kms 1 iBriigge n et al.|2005 



Heinz et al. 



2010). In addi- 



tion, sloshing of gas within the potential well of a cluster can give 
rise to similar velocities (Ascasibar & Markevitch 2006). Cosmo- 
logically generated turbulence is expected to be a relatively small 5 
to 15 per cent fraction of the thermal energy density in the cores of 
galaxy clusters ( |Lau et al.|2009l|Vazza et al.|2009||20iT) , 

The most direct method for measuring or placing limits on 
velocities on clusters is examining X-ray emission line widths. 
This was first done by examining a deep X-ray grating spectrum 
of A 1835 taken using the RGS (Reflection Grating Spectrome- 
ter) instruments on XMM-Newton (Sanders et al. 2010b). We fol- 
lowed this up with a study a sample of bright emission line objects 
(Sanders et al. 2011 1. As the RGS speedometers are slitless spec- 
trometers, the extent of the source contributes to the line width. For 
our most conservative analysis we assumed that the objects were 
point sources, and therefore obtained upper limits on the velocity 
broadening. For A 1835 we measured a 90 per cent confidence up- 
per limit of 274kms~ 1 and less than 500 kms -1 for several ob- 
jects in our sample. In A 1835 the upper limit corresponds to the 



fraction of energy in turbulence of less than 13 per cent of the ther- 
mal energy density. 

Other objects have recently been examined. |Bulbul et al.| 
([20T2J analysed XMM observations of A 31 12. Their most conser- 



de Plaa et al 



(2012) 



vative upper limit was around 900 kms 
recently examined two elliptical galaxies. Using resonance scatter- 
ing to provide lower limits and line widths to find upper limits they 
obtain 320 < V t iu- b < 720 km s" 1 for NGC 5044 and 140 < V turb < 
540^5-' for NGC 5813. 

In landers et al.||2010b} we went on to account for spectral 
broadening due to the spatial extent of A 1835. Using an Chandra 
image of the cluster in t he Fe-L band we red uced our limit to ~ 
180kms _1 . In addition, [Bulbul et al.](|2012} used XMM images 



to account for the surface brightness of A 31 12, finding improved 
limits of ~ 200 km s . 

When analysing the sample of objects in Sanders et al.| ( |20iT| >, 
we also applied an experimental technique of modelling RGS spec- 
tra using spatially-resolved Chandra maps of the sources, compar- 
ing the real line widths with predicted ones. This method had a 
number of systematic issues: using projected quantities, not includ- 
ing the uncertainties in the maps, assuming Solar abundance ratios 
and the finite sizes of the bins used. It has also been pointed out by 
M. Markevitch that we did not subtract our modelled line widths in 
quadrature from the observed line widths. 

In this paper, we improve on our previous analyses to take ac- 
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count of the spatial broadening in a robust manner. Firstly, we fit 
the RGS spectra, modelling the spatial distribution of the source. 
For some selected sources, we go on to use emission line profiles 
along the dispersion direction, subtracting continuum emission, to 
include the effects of the spatial broadening. We find a number of 
strong upper limits, but for a few sources we are able to find mod- 
erate detections of line broadening due to gas motions. 

The Solar abundance ratios of Anders & Grevesse ( 1989) are 
used in this paper. Velocity values quoted are the Gaussian width, 
a. 



2 SPECTRAL ANALYSIS 

The line width measured in an RGS spectrum includes the intrin- 
sic line width from random or turbulent velocities and thermal mo- 
tions. The thermal motion is included within the spectral model and 
is therefore accounted for during spectral fitting. The non-intrinsic 
contribution includes the instrumental line width (accounted for 
within the response) and the spectral broadening due to the extent 
of the source. The RGS instruments do not have slits, so the size of 
the source broadens the spectrum by 

AA^— AO A, (1) 
m 

where m is the spectral order and AO is the half energy width of 
the source in arcmin (note that the incorrect factor 0. 124 is given in 
some works, including Brinkman et al. 1998 ). In our previous anal- 
yses we placed an upper limit on the line width, putting a constraint 
on the velocity broadening. 

The contribution to broadening from the spatial component is 
a fixed broadening in wavelength. Velocity broadening gives a frac- 
tional broadening in energy or frequency. In principle, if we include 
a spectral model which characterises the fixed SX broadening from 
the extent, we can also determine the velocity broadening indepen- 
dently. If the data quality is too low we can constrain a combination 
of the size and velocity. 

In order to test whether this was possible we examined sim- 
ulated spectra. 1st and 2nd order XMM RGS1 and 2 response ma- 
trices were created with the SAS RGSRMFGEN tool. We supplied a 
custom Gaussian angular distribution with width a = 0.3 arcmin. 
Spectra were simulated with these responses in XSPEC (Arnaud| 
|T^96l > using a 1.5 keV BVAPEC model, 0.5 Z metallicity and 
including 500kms~' Doppler broadening. We took these simu- 
lated spectra and fit them using a point-source response matrix, a 
BVAPEC model, but convolved with a Gaussian with a fixed width 
in wavelength (parameterised by an angular scale, as in equation 
[Tp. We were able to recover the broadening velocity and Gaussian 
width of the source, cr, when fitting. 

2.1 Simple spatial modelling 

The data were prepared using the methods described in Sanders 
|et al.| ( |2011) . We used the same 90 per cent of the point spread 
function (PSF) extraction region in the cross-dispersion direction 
and 90 per cent of the pulse height distribution. This selection re- 
gion corresponds roughly to a 50 arcsec strip across the cluster in 
the cross-dispersion direction. Background spectra were extracted 
outside 98 per cent of the PSF. There have been some changes to the 
datasets included. We list the datasets used in Table [T] These data 
include observations of new objects such as MACS J2229. 7-2755 
(see Fig.[T](. We have removed some of the poorer quality datasets 



MACS J2229.7-2755 " 




Observed wavelength (A) 



Figure 1. RGS spectra of two new observations of MACS clusters in- 
cluded in this paper. The data where rebinned for MACS J2229. 7-2755 and 
MACS J0242.5-2132 to a minimum signal to noise ratio of 3 and 4, respec- 
tively. We note that MACS J0242. 5-2132 does not show the obvious narrow 
O VIII emission line as seen in MACS J2229. 7-2755. This appears to be due 
to a combination of a lower best-fitting oxygen metallicity (0. 15 ± 0.08 vs 
0.44±g;Jf Z ), broader emission lines and higher absorbing column density. 
The two clusters have similar temperatures, fluxes and redshifts. 

from the sample, in particular those where the background was a 
large fraction of the spectrum or the emission lines were weak. The 
datasets were reprocessed to take advantage of the updated cali- 
bration of the RGS line-spread function (XMM-CCF-REL-275). In 
addition, we increased the number of rows in the response matrices 
from the default value of 4000 to 16000. 

For most of the objects we fitted a single component BVAPEC 
spectral model (i.e. an APEC thermal model with variable abun- 
dances, thermal broadening and variable Gaussian broadening). 
The APEC model ( Smith e t al.|200T) we used was recreated from 
the APED 2.0.2 database with a temperature grid size of 0.01 dex, 
rather than using the default 0.1 dex tables. The PHABS photo- 
electric absorber component (Balucinska-Church & McCammon 
|1992[ l was applied to this model. We convolved the absorbed ther- 
mal model with a Gaussian of a constant wavelength width, param- 
eterised as the spatial scale. The data were fit in the XSPEC spectral 
analysis package, minimising the C-statistic during fitting. We fit 
the spectra between 7 and 28A for the first order data and 7 to 17A 
for the second order spectra. 

We note that using a spectral background extracted off-axis 
from the dataset itself may not be correct and depends on the source 
morphology and structure. If most of the flux from the object is 
emitted from within the extraction region, the background com- 
ponents will be correctly subtracted using off axis regions. If the 
source is more extended than the background region, the back- 
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Table 1. Details of the examined objects and datasets. Listed are the object names, redshifts, extraction position in decimal degrees (J2000), mean cleaned 
RGS 1 and 2 exposure and XMM datasets. 



Object 


Kedshirt 


RA 


Dec 


Exposure (ks) 


XMM observations 




2A 0335+096 


0.0349 


54.6713 


9.9669 


117.1 


0147800201 




A1068 


0.1375 


160.1855 


39.9529 


22.6 


0147630101 




A133 


0.0566 


15.6740 


-21.8805 


24.0 


0144310101 




A1795 


0.0625 


207.2182 


26.5937 


40.5 


0097820101 




A1835 


0.2523 


210.2580 


2.8787 


220.6 


0098010101 0551830101 0551830201 




A1991 


0.0587 


223.6318 


18.6449 


40.4 


0145020101 




A2029 


0.0773 


227.7339 


5.7446 


172.1 


0111270201 0551780201 0551780301 0551780401 












0551780501 




A2052 


0.0355 


229.1855 


7.0214 


116.7 


0109920101 0401520301 0401520501 0401520601 












0401520801 0401520901 0401521101 


0401521201 












0401521601 0401521701 




A2204 


0.1522 


248.1957 


5.5753 


85.9 


0112230301 0306490101 0306490201 
0306490401 


0306490301 


A2597 


0.0852 


351.3324 


-12.1243 


84.0 


0108460201 0147330101 




A262 


0.0163 


28.1924 


36.1528 


127.7 


0109980101 0504780101 0504780201 




A2626 


0.0553 


354.1272 


21.1467 


55.6 


00831502010148310101 




A3112 


0.075252 


49.4902 


-44.2381 


203.3 


0105660101 0603050101 0603050201 




A3581 


0.023 


211.8741 


-27.0183 


137.8 


0205990101 0504780301 0504780401 




A383 


0.1871 


42.0140 


-3.5291 


32.3 


0084230501 




A4059 


0.0475 


359.2542 


-34.7591 


46.6 


0109950101 0109950201 




A496 


0.0329 


68.4077 


-13.2616 


149.5 


0135120201 0506260301 0506260401 




AS1101 


0.058 


348.4948 


-42.7263 


98.0 


0147800101 




Centaurus 


0.0114 


192.2039 


-41.3125 


157.7 


0046340101 0406200101 




E 1455+2232 


0.2578 


224.3129 


22.3424 


36.1 


0108670201 




HCG62 


0.0137 


193.2738 


-9.2042 


156.9 


0112270701 0504780501 0504780601 




Hercules A 


0.154 


252.7841 


4.9924 


119.8 


0401730101 0401730201 0401730301 




Hydra A 


0.0539 


139.5249 


-12.0955 


93.5 


0504260101 




MACS J0242.5-2 132 


0.314 


40.6497 


-21.5406 


69.0 


0673830101 




MACS J2229.7-2755 


0.324 


337.4385 


-27.9270 


46.8 


0651240201 




MKW3s 


0.045 


230.4662 


7.7088 


32.6 


0109930101 




MS 0735.6+7421 


0.216 


115.4344 


74.2440 


55.2 


0303950101 




MS 2137.3-2353 


0.313 


325.0632 


-23.6612 


63.9 


0673830201 




NGC 1316 


0.00587 


50.6738 


-37.2079 


159.4 


0302780101 0502070201 




NGC 1365 


0.005457 


53.4016 


-36.1416 


111.6 


0505140401 




NGC 1399 


0.0046 


54.6212 


-35.4506 


126.0 


0400620101 




NGC 1404 


0.00649 


54.7156 


-35.5946 


41.9 


0304940101 




NGC 1550 


0.012389 


64.9080 


2.4100 


25.3 


0152150101 




NGC 2300 


0.007 


113.0821 


85.7093 


52.5 


0022340201 




NGC 3411 


0.0153 


162.6087 


-12.8451 


24.1 


0146510301 




NGC 4261 


0.00706 


184.8468 


5.8249 


94.6 


0502 1 201 01 




NGC 4325 


0.02571 


185.7776 


10.6210 


21.1 


0108860101 




NGC 4472 


0.00332 


187.4449 


8.0008 


83.6 


0200130101 




NGC 4636 


0.00313 


190.7080 


2.6876 


57.7 


0111190701 




NGC 4649 


0.00373 


190.9165 


11.5527 


76.0 


0502160101 




NGC 499 


0.01467 


20.7985 


33.4607 


49.1 


0501280101 




NGC 5044 


0.00928 


198.8498 


-16.3854 


133.7 


0037950101 0554680101 




NGC 533 


0.01739 


21.3808 


1.7592 


34.2 


0109860101 




NGC 5813 


0.0066 


225.2969 


1.7020 


162.3 


0302460101 0554680201 0554680301 




NGC 5846 


0.0057 


226.6223 


1.6049 


37.9 


0021540101 0021540501 




RBS 540 


0.0397 


66.4636 


-8.5601 


41.1 


0300210401 




RBS 797 


0.354 


146.8029 


76.3871 


17.7 


0502940301 




RXC J0605.8-3518 


0.141 


91.4749 


-35.3026 


24.1 


0201901001 




RXC J1044.5-0704 


0.1323 


161.1369 


-7.0691 


28.9 


0201901501 




RXC Jl 141.4-1216 


0.1195 


175.3518 


-12.2777 


29.5 


0201901601 




RXC J1504.1-0248 


0.2153 


226.0309 


-2.8044 


38.9 


0401040101 




RXC J2014.8-2430 


0.1612 


303.7154 


-24.5057 


26.4 


0201902201 




RXC J2149. 1-3041 


0.1179 


327.2822 


-30.7013 


26.8 


0201902601 




RXJ1720. 1+2638 


0.164 


260.0414 


26.6248 


69.2 


0500670201 0500670301 0500670401 




RXJ2129.6+0005 


0.235 


322.4164 


0.0886 


43.5 


0093030201 




Zw3146 


0.2906 


155.9152 


4.1865 


209.6 


0108670101 0605540201 0605540301 
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ground will contain source emission. The extraction region is a strip 
across the cluster and will contain emission from regions outside 
the core. Using an adjacent strip as background (as we obtain from 
the off-axis data) will remove much of this outer emission. This 
spectral component, as it comes from a more extended region will 
be more highly broadened than that from the source region. We are 
only interested in the best measurements from the centre, so it is 
favourable to remove the outer emission. This will particularly the 
case for a cool core embedded in a hotter medium. However, there 
are cases where the geometry of the source may make this subtrac- 
tion invalid, particularly if the source is not symmetric in the extrac- 
tion regions. For some objects it would be better to use background 
spectra from blank sky observations. We also redo the fits using 
template backgrounds generated with the RGSBKGMODEL tool to 
check the robustness of our backgrounds. 

We allowed the temperature, absorption, redshift, velocity 
broadening, spatial broadening, normalisation and some metallic - 
ities to be free in the spectral fits. The metals which we allowed 
to be free depended on the temperature of the object and whether 
there was any continuum in the spectrum. For some objects with 
high quality data or when required we fitted more than one broad- 
ened thermal component, to allow for more than one temperature or 
spatial scale. The parameters which were allowed to be free during 
the spectral fitting for each object are shown in Table|2] 

We examined the uncertainties in the parameters to the fits 
using both a conventional approach of varying the parameter un- 
til the fit statistic increases from the minimum by the required 
amount (referred to here as 8C) and using Markov Chain Monte 
Carlo (MCMC). Rather than use the standard Metropolis-Hastings 
MCM C algorithm, we us ed an affine-invariant MCM C sampler 
|Goodman & Weare|2010| as implemented in EMCEl£](Foreman- 
Mack ey et al.||2012) . We created a package based on EMCEE to 
analyse X-ray spectral models and data in parallel invocations of 
XSPEcj^] The advantage of an affine-invariant sampler over the 
Metropolis-Hastings algorithm is that it can sample complex pa- 
rameter space much more efficiently without any tuning of a pro- 
posal distribution. It uses a number of 'walkers' which simultane- 
ously examine parameter space within each step on the chain. The 
proposal distribution for a particular walker depends on the position 
of the other walkers. 

We started the walkers in a small cluster around the best fitting 
parameter values. 200 walkers were typically used, with a chain 
length of 2000 (per walker) following a burn-in period of 400. Re- 
peat fractions of the chains after burn-in ranged between 65 and 
77 per cent. For some objects, we found discrepancies between the 
MCMC and 8C statistic approaches or there were problems in the 
convergence of the chain. In these cases we used 4000 walkers, but 
decreased the length of the chain to 1000 and the burn-in period to 
200 to finish the analysis in a reasonable time period (Table|3]lists 
the number of walkers). When analysing the posterior probability 
distributions we use all of the walkers. The mean autocorrelation 
time of the velocity parameter for the walkers are typically less than 
3 per cent of the chain length and at worst 6 per cent of the chain 
length. We also examined the chains by eye to look for problems in 
the convergence. 

In Tableware presented the broadening results for each of the 
objects. The first two columns show the lo" and 90 per cent uncer- 
tainties (or limits) on the line-of-sight velocity obtained by allow- 



Table 2. Spectral fitting parameters used when fitting spectra. Nj is the 
number of temperature components used when fitting the spectra and Ng 
is the number of spatial components. If two spatial and temperature com- 
ponents are listed, two spatial components with independent temperatures 
were fitted. The listed metallicities were allowed to vary independently. If 
Fe is not listed, it was frozen at the Solar value because of a lack of contin- 
uum. Unlisted metals were tied to the Fe Solar ratio, except for He which 
was frozen at the Solar ratio. In addition to these parameters, we allowed the 
absorbing column density, normalisation (for each temperature component 
or spatial component), redshift and velocity width to vary. 

Object N T Ns Metallicities 



2A 0335+096 

A1068 

A133 

A1795 

A1835 

A1991 

A2029 

A2052 

A2204 

A2597 

A262 

A2626 

A3112 

A3581 

A383 

A4059 

A496 

AS 1101 

Centaurus 

E 1455+2232 

HCG 62 

Hercules A 

Hydra A 

MACS J0242.5-2132 

MACS J2229.7-2755 

MKW3s 

MS 0735.6+7421 

MS 2137.3-2353 

NGC 1316 

NGC 1365 

NGC 1399 

NGC 1404 

NGC 1550 

NGC 2300 

NGC 3411 

NGC 4261 

NGC 4325 

NGC 4472 

NGC 4636 

NGC 4649 

NGC 499 

NGC 5044 

NGC 533 

NGC 5813 

NGC 5846 

RBS 540 

RBS 797 

RXC J0605.8-3518 
RXC J1044.5-0704 
RXC Jl 141.4-1216 
RXC J1504.1-0248 
RXC J2014.8-2430 
RXC J2149.1-3041 
RXJ1720.1+2638 
RX J2129.6+0005 
Zw3146 



1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Fe Ni 

2 O Ne Mg Si Fe Ni 
1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 

2 O Ne Mg Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Fe Ni 

1 O Ne Mg Fe Ni 

1 O Ne Mg Si Fe Ni 

2 O Ne Mg Fe Ni 

1 O Mg Si Fe 

2 N O Ne Mg Fe Ni 

1 O Ne Mg Si Fe Ni 

2 O Ne Mg Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Fe Ni 

1 O Ne Mg Fe Ni 

1 N O Ne Mg Fe Ni 

1 O Ne Mg Fe Ni 

1 O Ne Mg Fe Ni 

1 O Ne Mg Ni 

1 O Ne Mg Fe Ni 

1 O Ne Mg Fe Ni 

1 O Ne Mg Ni 

1 O Ne Mg Fe Ni 

1 O Ne Mg Fe Ni 

2 O Ne Mg Fe Ni 

1 O Ne Mg Fe Ni 

2 O Ne Mg Fe Ni 

1 O Ne Mg Ni 

2 N O Ne Mg Fe Ni 
1 O Ne Mg Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Fe 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Fe 

1 O Ne Mg Si Fe Ni 

1 O Ne Mg Si Fe 

1 O Ne Mg Si Fe 

1 O Ne Mg Si Fe Ni 

1 O Mg Si Fe Ni 

1 O Ne Mg Si Fe Ni 



' See http : / /danf m. ca/em cee/| 

2 Available at https : //github. com/j eremysanders/xspec_emcee 
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Table 3. Velocity broadening results from the change in fit statistic (SC), both with standard off-axis and template backgrounds, and from the MCMC analysis. 
We show lo" and 90 per cent confidence uncertainties (kms~' ). The number of walkers used in the MCMC analysis is listed. We also show the temperature 
from the MCMC analysis, with 1 a uncertainties. If more than one temperature component was used, this is the value of the higher component. A blank MCMC 
result indicates that the limit was greater than 5000 km s~ 1 . 



Object 


5C 


5C 


5C Template 


SC Template 


MCMC 


MCMC 


MCMC 


Temperature 




(Iff) 


(90 per cent) 


(Iff) 


(90 per cent) 


(la) 


(90 per cent) 


walkers 


(keV) 


2A 0335+096 


< 190 


<300 


<210 


<330 


<200 


< 300 


200 


1.58 ±0.02 


A1068 


< 490 


< 710 


< 350 


< 570 


<480 


< 760 


200 


~ ^„^n 43 
z - o;, -0.38 


A133 


<420 


<580 


490+!«J 


490+lf« 


<380 


<540 


200 


3 17+0.59 
-'• -0.42 


A1795 


<550 


< 800 


<500 


<750 


<580 


<880 


200 


7 -ic+0.41 
-'"' J -0.32 


A1835 


< 240 


< 310 


< 220 


< 290 


< 200 


< 300 


200 


4 CQ+U.O/ 

^••>°-049 


A1991 


<390 


<460 


< 380 


<490 


zlu -170 


<400 


200 


1.68 ±0.06 


A2029 


<410 


<480 


<470 


<530 


770+ 70 
z/u -230 


<420 


200 


7 07+044 

- > -° z -0.29 


A2052 


<550 


< 700 


40()+220 
wu -340 


400 ±370 


<440 


< 640 


200 


1.72 ±0.03 


A2204 


<440 


<520 


110+ HO 
-""-230 


<510 




<480 


200 


3.68+°;« 


A2597 


460+1* 


<750 


<580 


< 680 


550+ 


•'•'"-510 


200 


Z * VZ -0.21 


A2626 
A262 2 


<470 
<400 


<670 
<570 


<690 
<590 


< 900 
<680 


<440 
<400 


<660 
<560 


4000 
200 


349+O.88 
1>fl -0.04 


A3112 


+ z u_ 140 


^ zu -270 


440+100 
^* u -150 


440+1 60 
#4U -310 


400+11° 

^ Ju -130 


43 n +150 
+JU_270 


200 


3 , , +0.13 
■ 5 - 11 -0.12 


A3581 2 


<300 


<460 


<400 


<590 


<280 


<420 


200 
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ing the fit statistic to increase by 1.0 or 2.71, respectively. These 
are followed by the same results when template background spec- 
tra are used. The next two columns show the same results obtained 
from the MCMC analysis. It is unclear how to calculate an error 
bar or limit from MCMC when a parameter is close to a hard limit 
(Okms -1 broadening here), relative to its uncertainty. Normally 
one would take the median value of a parameter and its percentiles 
from the chain. If the posterior probability distribution is significant 
at a limit then using a percentile to compute uncertainties will give 
finite values, rather than a value consistent with that limit. As many 
of our velocities are limits, we therefore used a different technique 
examining the marginalised posterior probability distribution. This 
distribution was computed for the velocity broadening in 250 linear 
velocity bins up to 5000 km s -1 . We selected the most likely ve- 
locity (this is the quoted value when significant). We moved down 
to lower values of marginalised probability until the required in- 
tegrated probability was obtained. The error bar that we quote are 
the minimum and maximum values of velocity that bound this in- 
tegrated probability region. 

We plot our values and uncertainties from the standard and 
MCMC methods for each object in Fig. [2] It can be seen that for 
most of the targets that we obtain good agreement between the two 
methods. However, there are a few cases in which they do not show 
consistent results. These include RBS797, where examination of 
the MCMC chain shows a long tail on the velocity posterior dis- 
tribution to high velocities and a secondary minimum in the fit 
statistic at high velocities. MS 0735.6+7421 shows a similar tail 
to high velocities, but the limit we can obtain is tighter. The tail on 
E 1455+2232 and RXC J2014.8-2430 are also longer than expected 
from the SC analysis. The two methods agree for RXCJ1504.1- 
0248 at the 1(7 level. Some objects show rather poor constraints 
(e.g. NGC499). Using template backgrounds makes little differ- 
ence to the results for most objects, although the change can be 
~ 200kms~ 1 . The worst cases include NGC499 and MKW3s, 
where the discrepancies are 600 — 800 km s~ 1 . 

In the determination of the broadening of the emission lines 
the main uncertainty is the extent of the source. We show in Fig. [3] 
the MCMC posterior probability contours of the velocity and spa- 
tial broadening for four objects. It can be seen that there is, of 
course, some degeneracy between the velocity and spatial broad- 
ening. However, this is not complete, as can be seen in A 31 12 and 
NGC 1399. 

In Fig. [4] we plot the temperature of each object against the 
best-fitting velocity (1<7 uncertainties, determining the quantities 
using MCMC). 



2.2 Incorporation of Chandra imaging 

The most obvious improvement to our method of fitting for the spa- 
tial extent of the source, is to use an X-ray image of the object from 
Chandra or XMM to measure the spatial broadening which should 
be included when modelling the spectrum. The profile of the image 
along the RGS dispersion direction is examined to compute how 
the spectrum is broadened. Either a wide band or narrower Fe-L 
band image could be used. X-ray CCD images have poor spectral 
resolution, so truly narrow band images cannot be made. Previously 
we used the RGSXSRC model in XSPEC to account for the spatial 
broadening in A 1835 ( |Sanders et al.|2010b| >. |Burbul et al.|p012) 
used MOS images and RGSRMFSMOOTH to account for source ex- 
tent in A 31 12. 



2A 0335+096 
A1068 
A133 
A1795 
A1S35 
A1991 
A2029 
A2052 
A2204 
A2597 
A2626 
A262 2 
A3 112 
A3581 2 
A3 83 
A4059 
A496 
AS1101 
Centaurus" 
E 1455+2232 
HCG 62 2 
Hercules A 
Hydra A 2 
MACS J0242.5-2132 
MACS J2229.7-2755 
MKW3s 
MS 0735.6+7421 
MS 2137.3-2353 
NGC 1316 
NGC 1365 
NGC 1399 
NGC 1404 
NGC 1550 
NGC 2300 
NGC 3411 
NGC 4261 
NGC 4325 
NGC 4472 
NGC 4636 
NGC 4649 2 
NGC 499 
NGC 5044 2 
NGC 533 
NGC 5813 2 
NGC 5846 
RBS 540 
RBS 797 
RXJ1720. 1+2638 
RX J2 129.6+0005 
RXC J0605.8-3518 
RXC J 1044.5-0704 
RXC J1141. 4-1216 
RXC J1504. 1-0248 
RXC J2014.8-2430 
RXC J2149.1-3041 
Zw3146 



5C (68%) 
SC (90%) 
MCMC (68%) 
MCMC (90%) 



100 



Velocity (km s" ) 



1000 



Figure 2. Error bars showing the 1 o~ (68 per cent) and 90 per cent confi- 
dence regions for each object. Results using 5C statistic and MCMC poste- 
rior probabilities are shown for each object. Objects marked with 2 include 
two spatial and/or thermal components in the spectral fit. 
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Figure 3. Confidence contours of velocity versus spatial broadening from 
the MCMC analysis. Contours show the 68.3, 90 and 95 per cent posterior 
confidence probabilities. 
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Figure 4. RGS temperature plotted against line width (using lo~ MCMC 
uncertainties). We plot the upper temperature if more than one temperature 
component is fitted. The lower curve shows where the energy density in 
random motions is 20 per cent of the thermal energy density. The upper 
curve shows the sound speed. 



2.2.1 Examining A 31 12 

Including a surface brightness profile in the modelling will fail if 
this profile does not accurately follow the emission lines that the 
RGS spectral fits are sensitive to. It might be thought that choos- 
ing a narrow band around Fe-L would give a good characterization 
of the line surface brightness. However, our testing has shown that 
even narrow bands often include too much continuum for many 
cool core clusters at larger radii, leading to over-correction of the 
line widths. The continuum at larger radii creates large wings in the 
response matrix which control the spectral fits. We demonstrate this 
in Fig. [5] which shows Chandra surface brightness profiles along 
the dispersion angle for an RGS observation of A 31 12. The pro- 
files are shown in a band tuned to the peak of the Fe-L emission 
(band B in Fig. [6) and the sum of two bands either side (bands A 
and C). The profiles are made in 1 arcsec bins along the dispersion 
angle, selecting the cross dispersion angle between —0.45 and 0.45 
arcmin. They are scaled by the number of counts in each band in an 
annulus from 1.4 and 1.9 arcmin radius. 

The surface brightness profiles continue to decline smoothly 
to large radii in both the Fe-L and continuum-like bands. If the 
profiles in the bands are subtracted (Fig. [5] bottom panel), the 
continuum-subtracted Fe-L emission is seen to be more compact, 
without the long tails at large radii. The strong Fe-L emitting region 
corresponds roughly to where the temperature drops below 3 keV 
in the cluster. 

We used the Fe-L 'On' (B band) profile to create response 
matrices using the 'angdistset' option to the RGSRMFGEN response 
generation tool. This option broadens the response matrix to ac- 
count for the broadening due to the dispersion axis profile given. 
We input the profiles between dispersion angles between —4 and 4 
arcmin. This was done for each of the A 31 12 datasets 0105660101, 
0603050101 and 0603050201, as the dispersion angles were dif- 
ferent in each case. The spectra and response matrices were com- 
bined with RGSCOMBINE and fitted in XSPEC. In this case we ob- 
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Dispersion angle (arcmin) 

Figure 5. (Upper panel) Chandra surface brightness profile along the dis- 
persion direction for the XMM 0603050101 observation for A 31 12. The 
profiles are shown in the rest frame Fe-L (B-band, 'On') and in the spectral 
region either side (A and C, 'Off') and are normalized at angles of 1.37 to 
1.9 arcmin. (Lower panel) Continuum-subtracted Fe-L emission, computed 
from difference between the B-band Fe-L and A+B profiles. 
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Figure 6. Chandra bands used for making line images for RGS analysis. 
Model spectra at example temperatures and 0.5 Z s are plotted. The effect 
of redshift were included. 



tained a best fitting broadening of Okms and an upper limit of 
175 km s~ 1 (90 per cent confidence). The C-statistic for the best fit 
in this case is 4148 with 4064 degrees of freedom. Our upper limit 
is roughly consistent with those of Bulbul et al. (2012), who found 
~ 200^8-'. 

Next, we instead used the continuum-subtracted Fe-L profile 
(bottom panel of Fig. |5j to create broadened response matrices, 
taking the values between dispersion angles of —0.4 and 0.4 ar- 
cmin, where the signal appears significant. A velocity broadening 
of 411 ±68kms~ 1 (lcr) is obtained using this dispersion profile. 
If we discard the negative values after the continuum subtraction, 
we obtain a value of 386kms~' instead. Both of these values are 
consistent with the confidence contours from the pure RGS spec- 




Offset (arcsec) 

Figure 7. Continuum-subtracted profiles in the cross-dispersion directory 
for various lines from A3112. We also plot the Chandra Fe-L continuum- 
subtracted profile convolved with the RGS2 cross-dispersion PSF at p = 
0.055. For O VIII in the 1st panel we extract the line from the wavelength 
range 20.35-20.52A and use both (20.05-20.22, 20.52-20.69)A as back- 
ground. For the 2nd panel we use 12. 54-12. 7lA as foreground and 12.72- 
12.90A as background. For the 3rd panel we use 11.75-12.36 and (11.37- 
11.50, 1 1.62-1 1.75)A as foreground and background, respectively. For the 
bottom 4th panel, we use 11.32-11.49A as the foreground and 11.14 to 
1 1 .3 1 A as the background. 

tral fitting shown in Fig. [3] The C-statistic for the best fit with this 
response matrix is 4103. When fitting for the spatial broadening 
(Section|2~Tjl, we obtained the same fit statistic. 

2.2.2 Systematic errors 

There are a number of possible systematic errors to consider. 
Firstly, the line width we measure is a function of several different 
emission lines, each of which may have different spatial broadening 
(and intrinsic width). Chandra X-ray images cannot separate these 
out spectrally. Even when not accounting for spatial extent, the sig- 
nal we measure will be some sort of average signal. Secondly, the 
radius at which we normalise the continuum may affect the result 
(particularly if there still line emission near the normalisation ra- 
dius). The band used to calculate the continuum may not be clean 
enough, for example if it contains emission lines. Another poten- 
tial problem is that the profile we input into the response generator 
needs to be truncated at a particular maximum offset, preferably 
where the tail of the surface of the brightness is zero but not includ- 
ing too much noise. Finally, the choice of background spectrum 
may have some effect, as in the previous section. 

The Chandra-derived continuum-subtracted dispersion axis 
profile can be compared to those extracted in the cross-dispersion 
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Table 4. Results using Chandra imaging. The Chandra observations used to make the images are listed, as are the bands used to make the line image (Line) 
and continuum image (Cont), the radial range where the surface brightness are normalised (Norm.) and the maximum offset of the profile along the dispersion 
direction when convolved with the response matrix. The background used for the spectral fitting (Bg) is either from the observation (O) or from a template 
(T). The velocities shown are limits or measurements with lo" or 90 per cent uncertainties. 'Selected background regions were chosen within these radii. 
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direction from the RGS observations themselves. However, unless 
the cluster is symmetric they should not necessarily be the same. 
Fig. [7] shows the cross-dispersion surface brightness profiles in 
four different lines, subtracting the continuum from adjacent spec- 
tral regions. Also shown on each plot is the Chandra continuum- 
subtracted Fe-L dispersion axis profile convolved with the RGS2 
cross-dispersion PSF The plot shows that for the three strongest 
lines, there is reasonable agreement with the Chandra profile. For 
Fe XXIII, the line appears narrower than the Chandra profile. This 
is probably because it is emitted from the cooler material in the 
cluster, located in the very central regions. This line is relatively 
weak and is unlikely to contribute significantly to the line width 
signal when fitting the spectrum. 

We can check how the continuum-subtracted profiles from 
Chandra depend on the choice of spectral bands. Fig.[8]shows pro- 
files where different continuum and line bands were used, and dif- 
ferent radii for normalisation of the continuum and line bands. It 
can be seen that the continuum band makes little difference to the 
continuum-subtracted line profile. The line band has some effect, 



where the two softer line bands have stronger peaks in the central 
region than the B band profile. The normalisation radius has little 
effect when the bands are close in energy (B and A+C), but does 
affect the result if they are further away in energy (B and D). 

We examined the effect of these possible systematic prob- 
lems by repeating the A3 112 analysis but varying some of the 
choices we made. Using a template background instead of the 
observation-derived background increased the best-fitting broad- 
ening by 30kms~' in this object. The line width decreased by 
lOOkms -1 when the dispersion axis angular range was increased 
from —0.4 to 0.4 to — 1 to 1 arcmin, using a template background. 
Increasing the angular range still further to 1.6 arcmin increased 
the line width slightly by 20 km s~ 1 . If we instead increased the 
radii where the line and continuum profiles were normalised from 
1.4 — 1.9 to 2.8 — 3.4 arcmin, this decreased the velocity by only 
6 km s~ 1 . The difference is more dramatic if a continuum band is 
chosen which is substantially different in energy from the line emit- 
ting band. Using band D as a continuum band, the broadening de- 
creased from 370 to 280 km s~ 1 when going from the inner to outer 
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Figure 8. Comparison of continuum subtracted profiles. The top panel show 
the effect of using different continuum bands (or none). The centre panel 
shows the effect of using different foreground bands. The bottom panel 
shows the effect of normalising the foreground and background in two dif- 
ferent radial ranges. The profiles have been normalised to have to have a 
total signal of unity. They have been smoothed using a running average, 
taking the mean of each point with its two neighbours. The key shows the 
line and continuum bands as a subtraction. 



normalisation radii. The inner region value is close to that when us- 
ing the adjacent A and C energy bands as continuum. The effect of 
the choice of band containing the line emission can be examined. If 
we combine the A, B and C bands, using D as the continuum, there 
is negligible change to the line width than when just using B. 

These tests show that the broadening we obtain can change by 
around 160kms _1 depending on the choice of bands, background 
and radii. This is the likely systematic uncertainty on any results 
of broadening. The range of values we obtain for A 31 12 (180 — 
500kms _I ) is consistent with the confidence contours shown in 
Fig.[3] 



2.2.3 Other targets 

We used this technique to include the surface brightness profile 
for several targets (listed in Table [4](. Unfortunately it cannot be 
straightforwardly applied to every object as certain conditions must 
be able to be met. Firstly, the core of the centre of the object must 
have strong line emission and the outer parts must not. This is 
achieved in cool core clusters. It should be unambiguous where 
the cool line-emitting region starts and the continuum region ends. 
We used Chandra temperature maps to help identify where to nor- 
malise the continuum and line profiles. Secondly, if the object has 
line emission in different emission lines which are spatially re- 



solved (for example, due to a temperature gradient over a large cool 
core) then the subtraction process often does not work, as the emis- 
sion lines move through different bands. This problem is seen in 
the Centaurus, A 2052 and A 3581 clusters, where the subtracted 
profiles are zero at large radii, become negative and then become 
positive. Finally, the object should have good Chandra data. The 
small PSF of Chandra means that we do not have to deconvolve 
it from the dispersion direction profile and enables us to easily re- 
move point sources. A future analysis using XMM data should be 
possible if its PSF is accounted for. 

For some distant objects without deep Chandra observations 
the bands A, B and C are rather narrow and contain too few counts. 
For these objects we therefore use A, B and C as foreground bands 
and a high-energy D band as a continuum band. For cooler objects, 
such as elliptical galaxies, the cool core is at lower temperatures. 
The Fe-L peak moves into the A band, instead of the B band. We 
therefore use the A and B bands as the foreground bands in these 
objects and band C as the background continuum band. The bands 
used for each observation are listed in Table [4] As with A 31 12, 
when convolving the continuum-subtracted emission line profile 
with the response matrices, we trim the profiles in dispersion angle 
where we do not see any signal above zero. These angular ranges 
are also listed in Table [4] 

Our velocity results are listed in Table|4] showing the statisti- 
cal uncertainties at the lcr and 90 per cent confidence levels. These 
measurements or limits were made using the SC approach. The 
same free parameters were used as when fitting in the previous sec- 
tion (except for the spatial broadening). 

There are a mixture of improved upper limits (such as A 1835 
and Zw3146) and possible detections of broadening of the emis- 
sion lines at 300 to 500 km s" 1 (e.g. A 3 1 12 and NGC 4636). 

We also investigate some of the possible systematic effects for 
the sample. For each object we compare the results using template 
and observation-derived spectral backgrounds. In some cases this 
can make large changes to the results (e.g. Hydra A, NGC 4636 and 
A 1068), but for the remainder the effect is small. For NGC 1399 
we investigate the choice of emission line band, normalisation ra- 
dius and extent of the profile included in the convolution. For this 
object the range in velocities is around 75 kms -1 , smaller than for 
A 31 12. We also examine the choice of continuum in A 2597 and 
the maximum dispersion angle included in the PSF convolution. 
Similarly to A 3 1 12, there is a range of ~ 150 km s~ 1 in the veloc- 
ity broadening, depending on the choice of parameters. 



2.3 Other methods 

The RGS instruments provide high spectral resolution images of 
the source in the cross-dispersion direction. The shape of the source 
in the cross-dispersion direction in a particular line can be directly 
measured from these images. If the source is close to symmetric 
in the dispersion and cross-dispersion directions, then the cross- 
dispersion profile can be used to account for the contribution of the 
extent to the line width. This method requires that there are enough 
photons in the lines for the extent to be measured (See e.g. Fig.[7jl. 
For many targets, this will not be the case for individual lines. 

A more complete approach would be to use Monte Carlo mod- 
elling of the cluster to fit all the observations, including imaging 
and the cross-dispersion information. Pete rson et al.| ( |2007) have 
used smoothed particles to model galaxy clusters, but have not yet 
used this technique to extract velocity width information. 
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3 DISCUSSION 

The measurements and limits on the velocity width we present are 
of the X-ray coolest gas, as measured using Fe-L and O VIII, char- 
acterised by gas at 3 keV or below. When ASTRO-H (Takahashi 
|et al.pOTO) is launched, its microcalorimeter will be most sensi- 
tive to the width of the Fe-K emission lines as it has fixed spectral 
resolution in terms of energy. Fe-K emission characterizes the hot- 
ter gas, which may give a wider range of velocities than the cooler 
material that the RGS is sensitive to. 

It is possible that the coolest gas may have sufficient inertia 
to avoid major velocity flows (or it would be destroyed by them). 
Previous results using RGS and Chandra have shown that the X- 
ray coolest gas consists of X-ray cool clumps embedded in a hotter 
medium ( Sanders et al. 2009, 20 10a I, providing that the two phases 
are in pressure balance. ASTRO-H will give complementary data on 
the hotter, volume-filling component. 

The turbulence or velocity structure of the coolest gas may be 
due to a number of different physical effects. The active galactic 
nucleus (AGN) may be generating kinetic feedback in its environ- 
ment, including jets and bubbling (Briiggen et al. 2005 
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20"T0]>. Cold fronts, if they require large flows (e.g. |Ascasibar & 



Markevitch 20061, may generate velocity broadening if they are 



within the spatial region examined by our analysis. We note that 
the physical region probed by our analysis is different for each of 
our objects. It is dependent on the redshift of the object as our ex- 
traction region is a fixed angular size and the physical extent of the 
line-emitted material within this aperture. 

The velocities we measure with RGS are comparable to, or 
exceed, the velocity width of optical emission line gas (after con- 
verted the quoted full width at half maximum, FWHM. to o~). Hatch 
eTaL] p007l > examined A 1068, A262, A496 and 2A 0335+096. 
Crawford & Fabian|(T992"l> ob served the objects AS 1 101, A 2597 
and A496. |Allen et al.| ( |1992| > measured line widths from A 1835, 
Zw 3146 and A 1068. 

Some of our results are inconsistent with results from previous 
analyses of line widths. Our measurement of 280 — 440 km s -1 for 
A 3 1 1 2 does not agree with the upper limit of 200 km s" 1 of|Bulbul 



|et al.| ( [2012| l when they included the effect of the spatial distribu- 
tion of the source. This is probably because the continuum was not 
subtracted from the surface brightness profile when computing the 
spatial broadening effect. In Sanders et al. (201 1 1 we used Chandra 
spatially resolved maps to model expected line widths, to subtract 
from observed line widths. This analysis appears to have overes- 
timated the spatial broadening effect in several objects. This can 
be seen by the disproportionate number of objects where the pre- 
dicted line width is wider than he observed line width, although 
some would be expected if the intrinsic velocity broadening is low. 
One likely contributing factor is the finite bin size of the input maps 
used to create the simulated spectra. In our previous paper the spa- 
tial and implied velocity broadening were also not added in quadra- 
ture to obtain the observed broadening. Our new method avoids the 
problems caused by our previous modelling using spectral maps. 

The velocities we obtain for NGC5044 and NGC5813 are 
consistent with those found by |de Plaa et al.|p012[ > using both line 
widths and resonance scattering. 

We note that for NGC 4636, the 430 - 580 km s _1 broadening 
we measure is apparently inconsistent with the 100 km s -1 upper 
limit inferred from the presence of resonance scattering (Werner 
|et al.| |2009:t. However, these values may still be in agreement if 
the broadening we measure is the result of ordered motion. The ra- 
dio source is likely to produce a bipolar outflow which will give a 



strong resonance scattering signal due to the low relative velocities 
of gas on either side. The emission lines, however, would be broad- 
ened by the difference in outflow velocities between each side. 
Alternatively, the disagreement may be because the spatial distri- 
bution of the line emission is not properly accounted for in our 
analysis. Although the Chandra surface brightness profile appears 
consistent with the location of the brightest and coolest part of the 
group (~ 0.6 keV temperature, in a 0.3 arcmin radius core, with 1 
arcmin wings), the profiles of individual lines may be different in 
this rather complex source. Another possible uncertainty is that the 
resonance scatter results may be affected theoretical uncertainties 
on the relative strengths of the Fe XVII emission lines as a function 
of temperature (see e.g. Sanders et al. 2008 1. 

In Section [2~Tj we measured the velocities assuming that the 
line surface brightness could be fitted well by a Gaussian model. 
There will be a systematic uncertainty associated with this assump- 
tion. In the case of A 31 12 this assumption is a reasonable one, 
as shown by the continuum-subtracted surface brightness profile in 
Fig. [5] and cross-dispersion line profiles in FigJT] The line width 
results for A 31 12 from this method (Table [3]and Fig. |3j also 
agree well with the results using the technique using the Chan- 
dra surface brightness profiles (Table R). Other objects with re- 
solved or marginally-resolved line widths show good agreement 
between the methods, including NGC 1399, MS 2137, A2204 and 
A 2597. The object with the largest difference between the methods 
is NGC 4636. This target has a complex morphology so the Gaus- 
sian assumption may not be good in this case. 

Detailed spatially-resolved line profiles (measuring line 
widths and positions as a function of radius) will enable jZhuravl-] 
|eva et al.|2012) the turbulent spectrum of the velocity field to be 
determined. We await the results from future X-ray missions flying 
high spectral resolution detectors, such as ASTRO-H, to examine 
line shapes in detail as a function of position. 

In our analysis, we used both a conventional and a MCMC 
method for examining uncertainties on parameters. For most of 
our objects, the results agreed. The use the affine-invariant MCMC 
sampler made it much easier to do the MCMC analysis compared 
to Metropolis-Hastings. With Metropolis-Hastings it is often diffi- 
cult to choose an appropriate proposal distribution, leading to the 
parameter space being poorly sampled because of too many or too 
few rejections. The algorithm also has difficulty when parameters 
are strongly correlated. The affine-invariant MCMC analysis was 
not a complete success, however, because in several cases we had 
to increase the number of walkers to get agreement between the 
conventional and MCMC results. Some parts of parameter space 
were missed in the MCMC analysis for some objects. Therefore, 
it was much easier to get quick, robust and stable results using 
conventional techniques rather than MCMC. The advantage of the 
MCMC approach was that we were able to examine the correla- 
tions of each parameter with every other and with the fit statistic 
(e.g. Fig. [3}. This allowed us to understand the modelling better. In 
a few cases, there are still divergent results between MCMC and 
the conventional technique, which may indicate that the C-statistic 
uncertainties are underestimated because of a tail in the posterior 
probability distribution. 



4 CONCLUSIONS 

We present limits or measurements of the velocity width of the 
soft X-ray emission for a sample of clusters and groups of galax- 
ies, and elliptical galaxies, using data from XMM RGS. In our 
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analysis we take two different approaches to including the spatial 
broadening intrinsic to the RGS instruments. For all targets, we 
model the spatial broadening with a Gaussian component in the 
spectral modelling. For a subset of objects we also use a second 
approach where Chandra imaging data in a continuum-subtracted 
line-emitting band is used to model the spatial broadening. 

Using Gaussian spatially modelling, for twenty of the tar- 
gets we find upper limits on the velocity broadening of less than 
500kms _1 and find five targets with limits of 300 km s _1 or lower. 
There are some targets for which we detect velocity broadening 
with 90 per cent confidence between 300 to 500 km s _1 , including 
A 31 12, HCG62 and NGC 1399. NGC4636 shows broadening if 
Chandra imaging is used in the spectral analysis. 



ACKNOWLEDGEMENTS 

ACF acknowledges the support of the Royal Society. We thank 
Randall Smith for supplying the APEC code. We are grateful to 
an anonymous referee for comments which significantly improved 
this paper. 



REFERENCES 

Allen S. W. et al., 1992, MNRAS, 259, 67 

Anders E., Grevesse N., 1989, Geochim. Cosmochim. Acta, 53, 197 
Arnaud K. A., 1996, in Jacoby G. H., Barnes J., ed, ASP Conf. Ser. 101: 

Astronomical Data Analysis Software and Systems V, p. 17 
Ascasibar Y, Markevitch M., 2006, ApJ, 650, 102 
Balucinska-Church M., McCammon D., 1992, ApJ, 400, 699 
Brinkman A. et al., 1998, in Proceedings of the First XMM Workshop on 

Science with XMM, http://xmm.esac.esa.int/ 
Briiggen M., Hoeft M., Ruszkowski M., 2005, ApJ, 628, 153 
Bulbul G. E., Smith R. K., Foster A., Cottam J., Loewenstein M., 

Mushotzky R., Shafer R., 2012, ApJ, 747, 32 
Crawford C. S., Fabian A. C, 1992, MNRAS, 259, 265 
de Plaa J., Zhuravleva I., Werner N., Kaastra J. S., Churazov E., Smith 

R. K., Raassen A. J. J., Grange Y. G., 2012, A&A, 539, A34 
Foreman-Mackey D., Hogg D. W., Lang D., Goodman J., 2012, 

arXiv: 1202.3665 

Goodman J., Weare J., 2010, Communications in Applied Math., 5, 65 
Hatch N. A., Crawford C. S., Fabian A. C, 2007, MNRAS, 380, 33 
Heinz S., Briiggen M., Morsony B., 2010, ApJ, 708, 462 
Lau E. T., Kravtsov A. V., Nagai D., 2009, ApJ, 705, 1 129 
Peterson J. R., Marshall P. J., Andersson K., 2007, ApJ, 655, 109 
Sanders J. S., Fabian A. C, Allen S. W., Morris R. G., Graham J., John- 
stone R. M., 2008, MNRAS, 385, 1186 
Sanders J. S., Fabian A. C, Frank K. A., Peterson J. R., Russell H. R., 

2010a, MNRAS, 402, 127 
Sanders J. S., Fabian A. C, Smith R. K., 201 1, MNRAS, 410, 1797 
Sanders J. S., Fabian A. C, Smith R. K., Peterson J. R., 2010b, MNRAS, 
402, LI 1 

Sanders J. S., Fabian A. C, Taylor G. B., 2009, MNRAS, 396, 1449 
Smith R. K., Brickhouse N. S., Liedahl D. A., Raymond J. C, 2001, ApJ, 
556, L91 

Takahashi T. et al., 2010, Proc. of SPIE, 7732, 77323J 
Vazza E, Brunetti G., Gheller C, Brunino R., Briiggen M., 2011, A&A, 
529, A17 

Vazza E, Brunetti G., Kritsuk A., Wagner R., Gheller C, Norman M., 

2009, A&A, 504, 33 
Werner N., Zhuravleva I., Churazov E., Simionescu A., Allen S. W., For- 

man W., Jones C, Kaastra J. S., 2009, MNRAS, 398, 23 
Zhuravleva I., Churazov E., Kravtsov A., Sunyaev R., 2012, MNRAS, 

422, 2712 



